Effective charge, collapse and the critical point of a polyelectrolyte chain 
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The charge of a polyelectrolyte (PE) controls myriads of phenomena in biology, biotechnology, 
and materials science, but still remains elusive from an understanding. Considering the adsorption 
of counterions on an isolated PE chain, an analytical expression for the effective charge - valid 
for all chain flexibility, for variable salt, in good solvents at all conditions, in poor solvents in the 
expanded (coil) state away from the critical point, and for gels, brushes and other PE systems in 
their expanded state - is derived. Phase boundaries and the critical point for the I-st order collapse 
transition, induced cooperatively by counterion adsorption and chain conformations, are calculated 
self-consistently. The size of the PE chain is found to be a single-valued function of charge. 
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The net charge of a polyelectrolyte (PE) chain in salty 
aqueous solutions is the most fundamental albeit elu- 
sive entity faced in systems containing charged polymers. 
Manifestations of charged polymers, whether m-RNA, 
ss- or ds-DNA, proteins, polysaccharides, or synthetic 
analogs such as polystyrene sulfonates, and the ways in 
which they fold conformationally and move in solutions 
depend crucially on their charge. These conformational 
fluctuations of PE chains under external stimuli are of 
critical importance to numerous biological and synthetic 
processes. A charged chain assumes in a good solvent 
a much larger dimension [j], Q than an uncharged self- 
avoiding-walk (SAW) chain due to electrostatic repulsion 
between like-charged monomers. This repulsive force re- 
sists the collapse of the PE in a poor solvent, but only 
up to a certain degree of hydrophobicity, beyond which 
it assumes similar to neutral chains a dimension way be- 
low the Gaussian value. [3, 4] Additionally, the PE has its 
counterions some of which adsorb (i.e., condense) on the 
monomers and modulate the effective charge. An opti- 
mization between electrostatic binding energy and trans- 
lational entropy controls this effective charge which, in 
turn, affects the size of the chain and vice versa. There- 
fore, a self-consistent procedure is necessitated to calcu- 
late the size and the effective charge of the PE molecule. 

A number of theories Q addressed the size dependency 
of a PE chain (or a PE gel@) on its effective charge, 
which is taken as a fixed parameter independent of its 
size. However, a fixed finite charge is incompatible to 
a collapsed PE chain in poor solvents as it leads to se- 
vere electrostatic energy penalty. Variable PE charge was 
addressed @ within the Oosawa model assuming rod- like 
chains. A three-state model[|| earlier self-consistently 
calculated the size and charge with the free counteri- 
ons partitioned into two domains - one within and the 
other outside the radius of gyration (R g ) of the PE. A 
successful analysis [9] of collapse due to counterion cor- 
relations in good solvents, however, failed to include the 
counterion adsorption effects and the right experimental 
trends. A theory that captures counterion adsorption, 
intra-chain electrostatic repulsion, and polymer confor- 
mational changes in a self-consistent way is yet to emerge. 



Experiments [T3, [H| anc ^ simulations [l^ have indis- 
putably demonstrated that a PE chain at a specific poor- 
ness of the solvent collapses to a globule. However, 
the stability of the exciting intermediate pearl-necklace 
phase[l]|, made of globules and elongated chain parts 
co-existing in a single chain, remains a controversial [l2j 
issue. In this letter we analyze the variational theory 
of counterion adsorption 0, [3], and derive an analytical 
expression for the charge of a flexible PE chain in good 
solvent conditions. Further, we determine the condition 
for collapse of the PE chain in a poor solvent and calcu- 
late the phase boundaries for this first-order coil-globule 
transition induced cooperatively by counterion adsorption 
and solvent poorness. Temperature controls both the sol- 
vent quality and the Coulomb strength. This formalism, 
fundamentally different from the Manning condensation 
theory [l5j, accounts for the interplay between chain con- 
formations and counterion adsorption from a neutraliz- 
ing plasma, as well as for the "oily" chemical nature of 
the polymer backbone. Addressing such interplay is cru- 
cial to understand the transformations between expanded 
and contracted polymer conformations encountered in a 
multitude of life processes, and to further explore a possi- 
ble decoupling of charge and size. The simplest candidate 
of coil-globule transition to exhibit the new paradigm is 
used to elucidate ideas pertinent to the swelling of ionic 
gels, behavior of charged polymers at surfaces and inter- 
faces, and assemblies such as viruses. 

The total free energy of the system consisting of an iso- 
lated flexible PE chain with ionizable monomers, coun- 
terions, and the added salt ions (all monovalent) in high 
dilution is taken to be[2| F = i*j where, 
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with the contributions being due to, respectively, the en- 
tropy of mobility along the chain backbone of adsorbed 
counterions (-Fi), the translational entropy of the free 
ions (F 2 ), fluctuations in densities of all mobile ions (F 3 , 
as in the Debye-Hiickel theory), adsorption (Coulomb) 
energy of the bound pairs (F4), and the free energy of 
the chain [T^j (F 5 ) resulting from its connectivity (con- 
formational entropy), the excluded volume (two-body), 
steric (three-body), and screened electrostatic interac- 
tions between monomers. Here f m = N c /N, a = M/N, 
with N, N c and M being the numbers, respectively, of 
monomers, of ionizable monomers, and of adsorbed coun- 
terions, ks the Boltzmann constant, and T the tem- 
perature. c s = cJq = UsIq/CI, and p = pl$ = NIq/CI 
with n s , Iq and fi being, respectively, the number of salt 
cations, the length of a monomer (Kuhn step length 
for a flexible polymer), and the volume of the system. 



K = y AitIb {p (fm — Oi) + 2c s } is the dimensionless in- 
verse Debye length. Here, R = kIq, and lg = C\/T is the 
dimensionless Bjerrum length with C\ = e 2 / AitkBtvdt), 
where e, eo, and e are the electron charge, vacuum permit- 
tivity, and the solvent dielectric constant, respectively. 
6 = (elo/eid) with e/, and d being, respectively, the 
local dielectric constant 0, [III, and the ion-pair sepa- 
ration. Further, (R 2 ) = Nlh = Nl 2 h = 6R 2 with, 

h - h/l and / = (JV C - M)/N = f m - a. o (a) = 
2/15(o <§; 1), l/3a(a ^> 1) is a cross-over function (a. H3| 
of a = k 2 jV!i/6. 

Generally, F = J^. Fi must be minimized self- 
consistently with respect to both / and l\ to obtain 
the equilibrium values of the variables. We note that 
/ features in F§ only in the fourth term that describes 
the long-range electrostatic interactions between charged 
monomers. However, the adsorption energy (F4) re- 
lated to the short-range ion-pair electrostatic attraction 
is more significant to F§ and shows a greater variation 
with / in most conditions. At low salt both the logarith- 
mic and the linear terms in F2 will dominate the / 3 / 2 
term in F3, and at high salt the variation of F3 with / 
is negligible. The full numerical analysis confirms these 
predictions, and allows us to propose an adiabatic approx- 
imation in which the chain free energy (F$) is decoupled 
from the rest. Then F a d = F\ + F2 + F4 are the relevant 
contributions that determine the effective charge /. This 




FIG. 1: Temperature induced collapse, phase boundaries for 
the chain charge and size, and the critical point of a fully ioniz- 
able polyelectrolyte chain. Fixed parameters are: N = 1000, 
p = 0.0005, 8 = 3.0, Ionizability f m = 1.0, 6 t = 400K, 
and c s i = 0.0. Critical values are: T* = 184. 9K, 103 = 8.5, 
ft = 2.3, /* = 0.020, and from T* , l* B = 4.867, w* = -1.163. 
Dotted curves are the analytical result [Eq. ©] in the ex- 
panded state. 



adiabatic approximation, however, progressively fails for 
lower molecular weights, higher Ib's and lower S's (F$ be- 
comes important and comparable to F4 in those cases), 
but, otherwise, remains applicable for a reasonably wide 
range of experimental conditions, especially at low and 
very high salts. Within adiabatic approximation, there- 
fore, .Fad does not contain the expansion factor and we 
can minimize over only a to obtain the effective charge. 
dF ad /da = gives us, using / = f m - a, 
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(fm = 1 for a fully ionizable chain). The present model is 
immediately related to experimental systems, where the 
"oily" backbone of the polymer presents heterogeneities 
in the dielectric function of the medium. Eq. © is the 
expression for the effective charge of an isolated PE chain, 
regardless of its flexibility, at all conditions in a good sol- 
vent, and in the expanded (coil) state away from the crit- 
ical point in a poor solvent. This closed form analytical 
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expression, as a function of the monomer density, maxi- 
mum ionizability, salt concentration, the Bjerrum length 
(temperature and the bulk dielectric constant), and the 
dielectric mismatch parameter S, is reported for the first 
time in literature. S — C2do has a non-universal constant 
C2 = (e/d)^ 1 , which can be determined from experiments 
using Eq. ^ . Further, the contributions to F ac j only de- 
pend on the number of ionizable monomers inside a poly- 
meric system, not on the size and shapes of the chains. In 
addition, the effective charge depends only on the aver- 
age monomer density p (in volume Q), but not explicitly 
on the molecular weight (N). The generality of the ef- 
fective charge expression warrants its validity for gels, 
brushes, and any other PE system in dilute solutions, as 
long as with changing / the variation in the adsorption 
energy for bound counterions (F4) is more significant to 
the variation of the intra-chain electrostatic interaction 
energy (term 4 in F§). For denser systems, higher coun- 
terion adsorption is expected in order to minimize the 
intra-chain repulsive energy. The size of the chain in the 
adiabatic (or expanded) state is determined by minimiz- 
ing F 5 with respect to h, using the / from Eq. ((2|). 

When the chain collapses (1% < 1, the Gaussian value) 
in a poor solvent away from the critical point, the chain 
collects most of its counterions reducing the effective 
charge and the electrostatic penalty, and the two- and 
three-body interaction terms dominate the free energy 
(F5). Here we employ the neutral chain approximation 
in which we assume / ~ or a = f m , and only F5 (with 
the fourth term in it being zero) survives. This approxi- 
mation breaks down (verified numerically) progressively 
closer to the critical point where the electrostatic interac- 
tions between monomers become comparable to the ex- 
cluded volume and steric interactions. However, away 
from the critical point dF$/dli = 0, in conjunction with 
l\ < 1 in the collapsed state, yields the expansion factor 
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This implies that R 2 g ~ R 2 g (9)k ~ Nl^N" 1 / 3 , or R g ~ 
TV 1 / 3 , which is the compact globule limit. Note that w 
must be negative to have a collapse (l± < 1). Now, ex- 
actly at collapse the binodal (equal chemical potential) 
states are, respectively, the / = (collapsed) and / 7^ 
(expanded). Note that for the present situation, chem- 
ical potential is the free energy per monomer. Equality 
of chemical potentials of the components in these two 
phases having zero osmostic pressure achieved in equilib- 
rium in a dilute solution yields, using Eq. j3| to deter- 
mine the free energy at collapse, 
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on T is given as w — (1 — 0t/T). Considering the tem- 
perature dependency of the Bjerrum length, Ib = C±/T, 
simultaneous solutions of T and a c in Eqs. |J2]) and (|4|) 
give us the temperature and the effective charge of the 
expanded chain at collapse. 
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where, a c satisfies Eq. @, and the dependency of w 



FIG. 2: Salt induced collapse of a fully ionizable polyelec- 
trolyte chain. Fixed parameters are: N = 1000, p = 0.0005, 
S = 3.0, Ionizability f m = 1.0, 9t = 400K, T = 300K, and 
therefore, Ib = 3.0, and w_= —0.333. Critical values are: 
5*! = 0.00305, w$ = 0.46, f* = 0.9, and /* = 0.051. Dot- 
ted curves are the analytical result [Eq. ^] in the expanded 
state. 

The above approximations which simplify the analysis 
by decoupling the chain free energy with rest and allow us 
to obtain the closed form expression of charge [Eq. ^] 
fail near the critical point where / is substantially higher 
than zero. Therefore, to obtain the critical point and the 
phase boundaries, we minimize the total free energy (sum 
of all five components) numerically with respect to both 
variables (i.e., / and The critical point and the phase 
boundaries for the temperature induced collapse of a fully 
ionizable chain in salt-free conditions are presented in 
Fig. [TJ The unique feature of this calculation is to take 
temperature as the only independent variable, thus ren- 
dering w and Ib as dependent functions through, respec- 
tively, the relations w = l- Q t /T, and l B = 18000/eT/ , 
assuming that the Bjerrum length in water at room tem- 
perature is ~ 0.75 nm, and the Kuhn length (Iq) is 
~ 0.25 nm [keeping in mind a polymer of type sodium 
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polystyrene sulfonate (NaPSS)]. The figures are represen- 
tations of phase transitions with the theta-temperature 
and dielectric constant as reference points, and the tem- 
perature a general experimental variable. Notable is the 
qualitative similarity of the charge (/) and size phase 
boundaries. There can not be a coexistence of expanded 
and collapsed states for the same effective charge of the 
chain. In other words, at equilibrium the chain can have 
only one size for one value of its effective charge. Fur- 
ther, the critical size of a fully ionizable chain at salt free 
conditions is larger than its Gaussian value (li = 1) due 
to intra-chain electrostatic repulsion, whereas for a neu- 
tral chain it is smaller (fx = 0.45 < f). Gradual reduc- 
tion of the ionizability (/ m ) leads to an increasing critical 
temperature and a decreasing critical size, for which the 
neutral chain values [H are recovered for f m = 0. The 
effect of ion- pairs @ that reduces the two-body interac- 
tion parameter only slightly without significantly affect- 
ing the phase boundaries is ignored. The phase boundary 
clearly suggests that counterion condensation and chain 
collapse are mutually cooperative processes in the case 
of a varying Coulomb strength. Lowering temperature 
in the expanded state helps counterion adsorption which 
reduces the effective charge and hence the resistance to 
collapse. Once the collapse occurs, charge becomes negli- 
gible with rampant counterion adsorption minimizing the 
electrostatic energy penalty. One may contrast this coop- 
erative behavior with a purely collapse induced collection 
of counterions [llj], where only the solvent quality (for the 
mixture of 1-propanol and 2-pentanone containing quar- 
ternized poly-2-vinylpyridine) is reduced keeping the the 



Coulomb strength (Bjerrum length times dielectric con- 
stant) virtually unchanged. Similar measurements of size 
and ionic conductivity as the temperature is varied would 
validate the major concepts presented in this letter. 

In Fig. [5] we present the phase boundaries for the salt 
induced collapse of a fully ionizable PE at 300K, at which 
the solvent is poor (note, 8t = 400!^). This shows that 
room temperature collapse of polymer of type NaPSS 
in water is possible in the presence of salt but unlikely 
at salt-free conditions (Fig. [T]) unless solvent poorness 
is drastically enhanced |ll|. There are remarkable quali- 
tative similarities between salt and temperature induced 
collapse (Fig.[TJ. The main ideas of the phase transitions, 
arising from the coupling between counterions and poly- 
mer conformations, are also valid for the volume transi- 
tions of PE gels and brushes. New experiments monitor- 
ing the counterion density inside the gels and brushes will 
stimulate a fundamental understanding of this important 
class of soft materials. One finally notes that computer 
simulations ignore the dielectric heterogeneity present 
in real experimental polymeric systems. Equivalently, a 
low value of 5 in our theory for which the counterion ad- 
sorption is weak raises the possibility of a necklace-like 
structure. Experiments 111], however, see rough spheri- 
cal objects with substantial incorporation of counterions 
within them, as predicted by our model. 
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